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Abstract 

We explore the role of soft gluon exchange in heavy quarkonium production at large 
transverse momentum. We find uncanceled infrared poles at NNLO that are not matched 
by conventional NRQCD matrix elements. We show, however, that gauge invariance and 
factorization require that conventional NRQCD production operators be modified to include 
nonabelian phases or Wilson lines. With appropriately modified operators, factorization is 
restored at NNLO. We also argue that, in the absence of special cancellations, infrared poles 
at yet higher orders may require the inclusion of additional nonlocal operators, not present 
in the NRQCD expansion in relative veloctiy. 

PACS numbers: 12.38.Bx, 12.39. St, 13.87.Fh,14.40Gx 



1 Introduction 

The production of heavy quarkonium offers a unique perspective into the process of hadronization, 
because the creation of the relevant valence partons, the heavy quarks, is essentially perturbative. 
Quarkonium production and decay have been the subject of a vast theoretical literature and of 
intensive experimental study, in which the effective field theory Nonrelativistic QCD (NRQCD) 
has played a guiding role. NRQCD offers a systematic formalism to separate dynamics at the 
perturbative mass scale of the heavy quarks from nonperturbative dynamics, through an expan- 
sion in relative velocity within the pair forming the bound state. In NRQCD, the description 
of the relevant nonperturbative dynamics is reduced to the determination of a limited number 
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of QCD matrix elements, accessible from experiment and, in principle, lattice computation. A 
characteristic feature of the application of NRQCD to production processes is the indispensible 
role of color octet matrix elements, which describe the nonperturbative transition of quark pairs 
in adjoint representation into quarkonia through soft gluon emission. 

An early success of NRQCD was to provide a framework for the striking Tevatron Run I 
data on high-p^- heavy quarkonium production [2], and it has been extensively applied to heavy 
quarkonia in both collider and fixed target experiments. A wide-ranging review of theory and 
experiment for quarkonium production and decay has been given recently in Ref. 0. Much of the 
analysis has been based on a factorization formalism proposed in jl], which offers a systematic 
procedure for the application of NRQCD to quarkonium production. It is fair to say, however, 
that, in contrast to quarkonium decay, fully convincing arguments have not yet been given for 
NRQCD factorization as applied to high-p^ production processes [SI E]. This omission may 
or may not be related to the current lack of confirmation for its predictions on quarkonium 
polarization at high px [S| • 

In this paper, we summarize progress toward the derivation of an appropriate factorization 
formalism for high-p^r quarkonia, illustrating our considerations with results on infrared emission 
at next-to- next-to-leading order (NNLO). At NNLO we find infrared divergences that do not fall 
precisely into the pattern suggested in Ref. p. These divergences may, however, be incorporated 
into color octet matrix elements by a technical redefinition that makes the latter gauge invariant. 
It is not clear whether this pattern extends beyond NNLO, and we conclude that NRQCD 
factorization must be examined further for production processes. In any case, all our results are 
consistent with the factorization of evolution logarithms in the ratio of momentum transfer to 
quark mass from nonperturbative matrix elements P]. 

In the results presented below, the relevant infrared divergence is proportional to f ^, where v is 
the relative velocity of the heavy pair in the quarkonium rest frame. The rotational invariance of 
this result (in the quarkonium rest frame) makes it possible to match the long-distance behavior 
of an arbitrary cross section to an octet matrix element in a manner that does not depend on the 
directions of energetic final-state gluons. In other words, we may factorize the perturbative long- 
distance contributions from the short-distance cross section, and replace them with a universal 
nonperturbative matrix element that has the same perturbative long-distance behavior, just as 
proposed in pn and extended in |n|. We begin our discussion with a brief review of NRQCD 
factorization at high transverse momentum. 

2 NRQCD in Fragmentation 

We discuss for definiteness the production of the J/ip and related heavy quarkonium states H 
in leptonic or hadronic collisions, A + B ^ H{pt) + X. To leading power in rric/pT, which we 
assume to be a small parameter, production proceeds through gluon fragmentation. According 
to conventional factorization theorems [7 , we have (keeping only the gluon) 

daA+B^H+xiPr) = daA+B^g+xiPr/zfi) ® DH/g{z,mc, fi) + 0{ml/p%) , (1) 
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where generally we pick the factorization scale to be of the order of px- In this expression, the 
convolution in the momentum fraction z is denoted by ®, and we have absorbed all information 
on the initial state into daA+B^g+x- If we also assume NRQCD factorization, we have in addition 
to (HD, 

dfTA+B^H+xiPr) = J2 dcrA+B^cc[n]+xiPT) (Cf ) , (2) 

n 

where the are NRQCD operators, classified by powers of relative velocity and characterized 
by the various rotational and color transformation properties of the cc state [n]. Assuming both 
^ and © to hold, we conclude that the gluon fragmentation function is related to the NRQCD 
matrix element by jH] 

DH/g{z, TTlc, /i) = X! dg^ccln]{z, /i, (Of ) , (3) 

n 

where dg^cc[n]{z, /i, ttLc) describes the evolution of an off-shell gluon into a quark pair in state [n], 
including logarithms of n/rric- 

In the following, we will study the fragmentation function itself, concentrating on infrared 
divergences at NNLO. First, however, we make some observations concerning the gauge trans- 
formation properties of NRQCD color octet matrix elements. 



3 Gauge Invariance and Wilson Lines 

Production operators for state H were introduced in Ref. jT] in the form 

O^{0) = xX^(O) {a^j^an) ^^/CU(0) , (4) 

where is the creation operator for state H, and where /C„ and /C^ involve products of color and 
spin matrices, and at higher dimensions of covaraint derivatives. Although the heavy (anti) quark 
fields (x) ip ^'I'e all at the same space-time point (here 0), the operator is not truly local, 
because the operator an creates particle H for out states, in the far future. In particular, 
operator- valued gauge transformations do not commute with the product a^f^an in general. 

A consequence of nonlocality is that the right-hand side of Eq. ^ is not gauge invariant 
in perturbation theory unless the individual product ip^]C'nX{^) and x^^n'^(O) are separately 
invariant. This is the case when the /C^'s specify color singlets, but not when they specify color 
octets. A related issue arises in the field-theoretic definitions of fragmentation functions P], such 
as Dff/g in Eq. (jT)) above. It is resolved by supplementing the bi-local products of fields by 
nonabelian phase operators, or Wilson lines: A] = exp [—ig /o°° dX I ■ A{x + XI)]. In contrast 
to parton distributions, moments of fragmentation functions do not result in expectation values of 
local operators [IQl, and the corresponding Wilson lines are not guaranteed to cancel identically. 
For a color octet combination, the gauge field is given in adjoint representation, just as for a 
gluon fragmentation function, for which the operator ^i{x)ab multiplies the gauge covariant field 
strength Fj^'^ . To be definite, and for convenience in relating NRQCD operators to the gluon 
fragmentation function, we will choose / as a lightlike vector in the minus directon: = 
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Our gauge invariant redefinition of production operators in octet representation is now found 
by tlie replacement (which we refer to as a gauge completion), 

O^{0)^X^}Cn,cm^][0,AU[aUH) H'^,A],aX^K,am, (5) 

where we have exhibited the color indices of the octets. In perturbation theory, the Wilson 
lines generate propagators of the form ■ k + ie) when they carry gluon momentum k, and 
the gluons couple to the Wilson line at vertices Qsl'^Cahc, with Cabc structure constants. These 
are precisely the same as the propagators and vertices of the eikonal approximation for the 
emission and absorption of soft gluons by an energetic gluon of momentum . We will use this 
correspondence below. 

It is worth noting that the replacement (0) isn't really necessary if gauge dependence in 
the matrix elements and coefficient functions is infrared safe. This appears to be an implicit 
assumption in the discussion of Ref. jT]. As we shall see, however, although gauge dependence 
is certainly free of infrared singularities at next-to-leading order, at NNLO this is no longer the 
case. 

4 Eikonal Approximation for g ^ H 

For this discussion, we study the perturbative expansion of the fragmentation function for g 
cc[no\ + X with cc[?t,o] = c(P/2 + q)c{P/2 — q), a pair with total momentum P and relative 
momentum q, always projected onto a color singlet state. We can, of course, further project 
various spin and orbital angular momentum states [n] as in Eq. (jS)). 

Beyond lowest order, the computation of these diagrams tests NRQCD factorization, Eq. Q, 
which is verified to the extent that we can absorb, that is match, any infrared divergences into 
the perturbative expansion of NRQCD matrix elements 11 . We will be interested in gluons 
whose energies are much below rUc- For NRQCD factorization to hold these soft gluons must 
factorize from the higher energy radiation that generates logarithms of pt/ttLc. ^ 

Because we need only the infrared structure of the diagrams from gluons of very low momen- 
tum, we may suppress overall color and combinatorial factors, as well as momentum factors that 
depend only on the scale of m^, including the LO fragmentation of the parent gluon, which is 
off-shell by /i^ > 4m^, into the quark pair. 

Figure la shows the lowest order diagram g cc, and Figs. lb,c show typical contributions 
with single soft gluon emission in the fragmentation function, all in cut diagram notation. Both 
Fig. lb and Ic represent a class of four diagrams, in which gluons are attached to one of the two 
lines of the quark pair in the amplitude and its complex conjugate. In each diagram we may 
treat the quarks in eikonal approximation, in effect replacing them with time-like Wilson lines, 
whose propagators and vertices are given hy i/ {(3 -k + ie) and —iggTaP^, respectively, for quarks, 
with Ta a color generator in fundamental representation and (3'^ the time-like quark four- velocity. 
The corresponding approximation for antiquarks simply changes the sign at the vertex. 

^We will show elsewhere that this is indeed the case. 
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(a) (b) (c) 

Figure 1: a) Lowest-order fragmentation function for g —>■ cc. There are no interactions on the 
eikonal quark pair or the Wilson hne that corresponds to an eikonal gluon of four-velocity /. b,c) 
Representative NLO contributions to g —>■ cc fragmentation in eikonal approximation. In these 
figures the parent gluon is contracted to a point, represented by the dark circle, because it is 
off-shell by order rric 



The vertical lines in Fig. lb and Ic represent the quark pair final state, including a pro- 
jection onto a color singlet configuration. The full fragmentation function is found by cutting 
the remaining lines in all possible ways. When the gluon is cut in Fig. lb, in particular, the 
resulting diagrams have an uncancelled infrared pole in dimensional regularization, beginning at 
order v"^ ~ q'^/m'^ when the gluon momentum k ^ 0. This gluonic contribution, however, which 
involves only interactions between the pair, may be matched straightforwardly to a color octet 
matrix element, with which it shares an identical topology [HIE]- 

The divergent part of Fig. lb may be isolated by evaluating the four relevant diagrams in 
eikonal approximation and then expanding them in f = q/rric- Alternately, we may expand 
the momentum-space integrand, isolating the f ^ term before the integration. In D = A — 26 
dimensions, and in the quark pair rest frame, the latter approach leads readily to an expression 
in terms of the total and relative momenta P (P^ = m^) and q, respectively, 

S(i^)(P,g) = IQg^^^^^^ j ^-^^5{e)[q,{P-k)-{q-k)Pj\ 

X [q^iP-k) - iq-k)P''] ^ 



212 



[(P ■ k) 

16 as q^ 1 , . 

= yvp^- + ---' 

where we have exhibited the infrared pole, regularized for e < 0. This infrared pole, along with 
an appropriate color trace, appears at NLO as a multiplicative factor times the lowest-order 
fragmentation function, as found, for example, in Refs. [^1 112j. 

In the integral, the terms in square brackets correspond to the lowest-order vertices for the 
operator q'^F^pP^ ■, with F the field strength. In the quarkonium rest frame, this corresponds to 
an electric dipole transition from octet to singlet color states [H]. An advantage of this expansion 
is that, because the field strength decouples from scalar polarized lines, all singularities associated 
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with gluon momenta collinear to l^^ vanish on a diagram-by-diagram basis, and double poles in e 
are entirely absent in the calculation. By the same token, the integral in (jH)) is gauge invariant 
in the polarization sum for the final-state gluon of Fig. lb. 

Finally, in Fig, Ic, the soft gluon connects the pair to the Wilson line. In covariant gauges, 
this diagram has double and single infrared/collinear poles in dimensional regularization. These 
poles, however, cancel against analogous contributions from the corresponding virtual diagram. 
Imaginary parts of the virtual diagrams, of course, cancel in the fragmentation function, which 
is real. In practical terms, the contribution of Fig. Ic plus its virtual counterpart contributes 
solely to the perturbative factor dg^cc[n]{z, fi,mc) in Eq. Q. 

We note that the infrared behavior of all of the diagrams of Fig. 1 is common between the 
fragmentation function and a generic cross section in which a high-p^ gluon of momentum / 
recoils against the heavy quark pair. Indeed, the same cancellation mechanism for the infrared 
divergences of Fig. Ic is referred to specifically in Ref. PP as an essential step in the justification 
of NRQCD factorization for production cross sections. To the extent that all the infrared poles 
of diagrams that do not share the topology of NRQCD matrix elements cancel, the matching 
of perturbative infrared poles with those of the effective theory is ensured. The result of this 
cancellation is referred to as topological factorization [S]. The same considerations that lead 
to the inclusion of nonabelian phases to the fragmentation function, however, suggest that the 
cancellation of non-factored soft gluons may be nontrivial, and should be checked beyond NLO. 
Of course, a full NNLO calculation of either cross sections or fragmentation functions would be 
daunting. Fortunately, the analysis of the relevant infrared behavior at NNLO requires only the 
eikonal approximation, and is therefore a much more manageable, although still extensive, task, 
to which we now turn. 

5 Fragmentation Function at NNLO 

We can readily generalize our NLO discussion to NNLO. Once again, because we are interested 
in gluons whose energies are much less that rric, we can adopt the eikonal approximation, and 
neglect all dynamics at the scale of the quark masses. In NRQCD, these approximations are 
relevant to gluons of order rricV, with v the relative velocity, but of course in the perturbative 
calculation the gluon momentum goes to zero. We note that for octet-singlet transitions at 
NNLO, we need not include virtual gluon exchange between the quark and antiquark, except for 
diagrams that are already topologically factorized. As a result, we will not need not consider 
violations of the eikonal approximation from momentum regions at or below the scale rricv'^, or 
the resulting 1/v singular behavior. 

Representative diagrams for the fragmentation function are shown in Fig. 2. The full infrared 
fragmentation function is again generated by taking all allowed cuts of the remaining lines of 
each such diagram. 

We are concerned only with diagrams that connect octet to singlet quark states, and which 
are not topologically factorized, since these are the potential sources of nonfactoring behavior in 
both the fragmentation function and related cross sections. We recall once more that the familiar 
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argument for NRQCD factorization is based on the conjecture that all infrared regions in these 
diagrams cancel after this limited sum over cuts p. In fact, we shall see that this is the case 
at NNLO only if we employ the gauge-completed definitions for NRQCD matrix elements, as in 
Eq. (0) above. 




(a) 



(b) 



(c) 



Figure 2: Representative NNLO contributions to g ^ cc fragmentation in eikonal approximation. 

As at NLO, we will concern ourselves here with infrared behavior at order v"^, which may be 
computed by evaluating each diagram in eikonal approximation and then expanding their sum in 
q/P, or by expanding the diagrams first, as in Eq. Noting that the two approaches give the 
same answer, we summarize the results here, and provide details of the calculations elsewhere 

The individual classes of diagrams in Fig. 2a and 2b, for which two gluons are exchanged 
between the quarks and the Wilson line, satisfy the infrared cancellation conjecture of Ref. [T], 
by summing over the possible cuts and connections to quark and antiquark lines, as do diagrams 
that have three gluon-eikonal vertices on the quark pair and one on the Wilson line. For the 
class of diagrams related to Fig. 2c, however, with a three-gluon interaction, this cancellation 
fails. Expanding again to second order in the relative momentum q, the full contribution from 
Fig. 2c, found by cutting the gluon line ki and the Wilson line can be written by analogy to Eq. 
(la as 



S(2^)(P,g,/) 



16i g%^' 



277 5{kl) l^K^x[kuk2\ 



X [g^(P • k,) - {q ■ k,)P>^] [q'^iP ■ k,) - (g ■ k^)P-\ 
1 



X 



X 



[P ■ fci + ief [P-k2- ieY 
1 



[kl-ie] [{k^-k^Y-ie] [l-{k^-k2)-ie] 



(7) 



where V„^\[ki, ^2] represents the momentum part of the three-gluon coupling. As in Eq. we 
have suppressed color factors and momentum dependence at the scale rric. 

As observed above, the field-strength vertices eliminate coUinear singularities on a diagram- 
by-diagram basis. The leading singularities in (|7j) and related diagrams are therefore never worse 
than l/e^. Summing over all such contributions, however, we find a noncancelling real infrared 
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pole in the fragmentation function, which may be written in invariant form as 






(8) 



In the rest frame of heavy-quarkonium (P = 0), this becomes simply 



S(P,g,/) 



3e 4m2 




(9) 



where v is the relative velocity of the heavy quark pair. We will give the extension of this result 
to all powers in v elsewhere [T^ . 

Eq. shows explicitly the breakdown of the simplest topological factorization of infrared de- 
pendence at NNLO. Its presence implies that infrared poles would appear in coefficient functions 
at NNLO and beyond when the factorization is carried out with octet NRQCD matrix elements 
defined in the conventional manner, Eq. Q. On the other hand, when defined according to its 
gauge-completed form (jSj), each octet NRQCD matrix element itself generates precisely the same 
pole terms given in (jH)) above. Thus, at least at NNLO and to order f^, NRQCD factorization 
can accomodate these corrections. We note, however, that at this order the NNLO correction is 
independent of the direction I of the Wilson line, while Lorentz invariance alone would seem to 
allow a correction of the form (g-Z)^/P^. The presence of poles with this coefficient could indicate 
problems with factorization, that is, with matching the poles of the fragmentation function to 
those of the perturbative cross section. These could arise when the integration over q is not 
rotationally invariant, as might happen for polarized final states, and/or final states with gluons 
of energies of order rric in the rest frame of H, since their directions are arbitrary. We are unable 
at this stage to rule out the occurance of such terms at higher orders in the strong coupling, or 
in higher powers of v"^. 

6 Summary and Future Studies 

We have discussed the gauge invariance properties of color octet matrix elements in NRQCD 
as they appear in fragmentation functions. Gauge invariance and factorization require that 
they include nonabelian phase operators to match otherwise nonfactoring infrared divergences 
beginning at NNLO. 

We have shown that at NNLO there are infrared (not collinear) poles from fragmentation dia- 
grams that are not topologically equivalent to conventional NRQCD operators. We have observed 
that rotational invariance in the quarkonium rest frame is an important consistency condition for 
the possibility of NRQCD factorization in terms of the gauge-completion of NRQCD operators. 
Should rotational invariance fail at higher orders, we anticipate that it will be necessary to com- 
plement the NRQCD classification of nonperturbative parameters for quarkonium production 
with matrix elements of additional operators, involving nonlocal products of the field strengths 
encountered above. 
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Finally, although our explicit calculations are carried out in the eikonal approximation for 
gluons of energies far below m^, it is not difficult to verify that these gluons do not interfere with 
the generation of standard evolution for the g —>■ H fragmentation functions. The arguments for 
this result will be given elsewhere, along with details of the calculations described above. 
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